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We examine several effects of an applied magnetic field on Anderson-type models for both single- 
and two-level quantum dots, and make direct comparison between numerical renormalization group 
(NRG) calculations and recent conductance measurements. On the theoretical side the focus is on 
magnetization, single-particle dynamics and zero-bias conductance, with emphasis on the univer- 
sality arising in strongly correlated regimes; including a method to obtain the scaling behavior of 
field-induced Kondo resonance shifts over a very wide field range. NRG is also used to interpret 
recent experiments on spin-^ and spin-1 quantum dots in a magnetic field, which we argue do not 
wholly probe universal regimes of behavior; and the calculations are shown to yield good qualitative 
agreement with essentially all features seen in experiment. The results capture in particular the 
observed field-dependence of the Kondo conductance peak in a spin-i dot, with quantitative devi- 
ations from experiment occurring at fields in excess of ~ 5T, indicating the eventual inadequacy of 
using the equilibrium single-particle spectrum to calculate the conductance at finite bias. 
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I. INTRODUCTION 

Understanding electronic transport in quantum dots 
remains a major challenge for theorists working on cor- 
related electron systems. Conductance at low energies is 
often dominated by one of a number of Kondo effects,^ 
in which the strong localized interactions on the dot(s) 
induce non-trivial many-body physics. Over the years a 
wide range of such Kondo effects have been predicted and 
observed, in single and multiple quantum dot devices of 
various geometries 3"— 

Here we consider a single quantum dot, tunnel-coupled 
to two leads^r— in an effective one-channel fashion. 
While the dot will in general hold many electrons in its 
quantized levels, only those close to the Fermi level con- 
tribute in practice to electronic transport provided the 
mean level spacing is sufficiently large, and the rest can 
be neglected with relative impunity. Typically just one 
level is important, but occasionally one observes the case 
of two relevant levels, where the physics is richer; in- 
cluding e.g. a quantum phase transition between Fermi 
liquid and underscreened Kondo phase»i2r— which has 
been observed in several experimental guises j^^-22, 

We present and examine critically a number of results 
falling under the umbrella of magnetic field {B) effects in 
these single- and two-level quantum dots, the appropri- 
ate models for which are specified in sec. |lTl The paper 
consists of two related parts. In the first (sees. IIIIIIV|) . 
using mainly Wilson's numerical renormalization group 
(NRG) method)^"— we consider magnetization, single- 
particle dynamics and the zero-bias conductance, with 
emphasis on the universality and scaling behavior aris- 
ing in the strongly correlated regimes of the models. 

Even for single-level quantum dots described by an 
Anderson impurity model,— i^ there are still open ques- 
tions regarding single-particle dynamics in the presence 
of a magnetic field; our primary concern being the field- 



induced Kondo peak splitting in the equilibrium single- 
particle spectrum. This has been analyzed by a number 
of authors and techniques'^"— but the results are not in 
complete agreement j^^i^i^ld^ We show that there exists 
an algorithm by which NRG can obtain the universal be- 
havior over many orders of magnitude of field strength 
but that, eventually, even the most accurate NRG calcu- 
lations cannot completely resolve the universal splitting 
at very large fields. The corresponding situation for the 
two- level model is considered in sec. lIVB] We also obtain 
the field and temperature (T) dependence of the zero- 
bias conductance, and for T = in particular generalize 
the Luttinger integral analysis of ref. [23 to encompass 
a finite magnetic field; leading to an exact result for the 
conductance for any field, and insight into the rather sub- 
tle differences between the limits B = Q and S — > for 
the underscreened triplet phase of the model. 

In the second part of the paper (sec. |V| we turn to 
comparison with experiment. Two recent sets of conduc- 
tance measurements on quantum dots in a magnetic field 
are considered, from the groups of Kogar^ (on an ef- 
fective one-level dot) and Goldhaber-Gordon^ (on both 
effective one- and two- level systems). From comparison 
to NRG results, we are able to determine reliable bare 
model parameters for the Anderson-type (as opposed to 
Kondo) models considered in sees. |ll] - IIVI as relevant 
to experiment. With these, our NRG calculations are 
shown to yield very good qualitative agreement with es- 
sentially all features observed in both experiments4^i^ 
In particular, we show that theory can in fact explain the 
evolution of the field-induced splitting of the Kondo con- 
ductance peak observed in ref. |4J - including a simple 
explanation for an observed crossing in the peak split- 
tings of two different quantum dots. The agreement is 
essentially quantitative up to field strengths of around 
a couple of Kondo scales, but beyond that our calcula- 
tions deviate from the experimental data. This reinforces 



results from a recent study using the scattering states 
NRG^*^ and earlier renormalized perturbation theory and 
NRG calculations,—!^ showing that the commonly used 
approximation of calculating the source-drain bias de- 
pendence of the conductance from the equilibrium spec- 
trum is unsuitable for making quantitative comparisons 
to experiment sufficiently far out of equilibrium. Indeed, 
until more progress in non-equilibrium theory is made, 
we suggest that experiments should instead aim to make 
comparison with the magnetic-field dependence of the 
zero-bias conductance. 



II. MODELS 

Each model considered in this work consists of a single 
interacting quantum dot region, tunnel coupled to a pair 
of non-interacting metallic leads. As mentioned above, 
we focus on the situation where the mean level spacing 
of the dot is sufficiently large compared to the dot-lead 
tunneling strength that generally only one, or occasion- 
ally two, levels are involved in transport. 

When just one dot level is relevant the standard model 
is the Anderson impurity model (AIM)<ii2ii2^ Here the dot 
itself is described by 



^^™ = E (^1- + lUn^-^ 



(1) 



where nio- = d[^di^ counts the a spin electrons on the 
dot level, U is the on-level Coulomb replusion/charging 
energy, and eia = ei — crh the level energy. The latter 
includes a Zeeman coupling to an external magnetic field 
B with h = ig/iB^ and a = +/— for ^ / |-spin electrons. 
In the case of two relevant dot levels, the dot Hamiltonian 
is naturally more complex. We choose to work with the 
following two-level model (2LM) 



H. 



2LM 



i,(T 



2Uf>■^-a]n,„ 



U'n^n2 — Jh^i ■ Sj 



(2) 
which has previously been shown to capture the key 
physics of two-level quantum dots in the absence of a 
magnetic fieldi^ Here n^ = J^a ^ia i^ the total number 
operator for level i (=1,2), and s^ is the local spin oper- 
ator with components sf = J2a.a' Aa'^aa'^^icr' ("'crcr' are 
the Pauli spin-i matrices). In addition to the on- level 
Coulomb repulsion U (taken to be identical for levels 1 
and 2 for simplicity), the model includes an interlevel 
Coulomb repulsion U' plus a ferromagnetic (Hund's rule) 
exchange coupling of the spins of the two levels, param- 
eterized by Jh- 

In each case, the dot Hamiltonian is supplemented 
by coupling to two equivalent, noninteracting 'left' 
and 'right' leads, themselves described by Hi^ = 
El. Ek,^ ekaL^«k^<T (^ = ^'^)' where the most gen- 
eral tunnel coupling to the leads is of form Ht = 



J2u J2i,k,a ^iki^i4a^kua + ^-C.) (the sum over level index 
i involving just i = 1 in the case of the AIM). The L and 
R lead chemical potentials are /xz, and fj^R respectively, 
such that for fi^ ^ fJ-B. a non-zero current flows. 

Analyzing the interacting models out of equilibrium 
is a formidable task (see e.g. ref. |43 for a recent dis- 
cussion) and in practice we consider the equilibrium sit- 
uation. This has an immediate benefit, for the AIM 
Hamiltonian then reduces exactly to an effective one-lead 
model by defining c,^^ = iVikLai^j^„ + V.ii^Rai^j^^)/Vik with 
V^ = V^^ + V^j^ (with i = 1), since the corresponding 
orthogonal combination of lead states is entirely decou- 
pled from the dot. The two- level dot Hamiltonian un- 
der this transformation does not generally separate so 
pristinely: excepting the special case of Viki — ^ik cos 9, 
Vi\^R = Vi-k sin 6*, the dot remains coupled to two leads J^ 
However, over a wide range of parameter space the sec- 
ond lead couples sufficiently weakly that it may in prac- 
tice be neglected on energy scales of practical interest.— 
As such, for both the AIM and 2LM we work with the 
effective one-lead description embodied in 

N 

Hi^+Ht = E £kcLck.+E E ^.k(4ck.+H.c.). (3) 



k,cr 



i—1 k.cr 



with iV = 1 for the AIM and TV = 2 for the 2LM. 

We consider the standard casei of a symmetric, flat- 
band lead of half-width D and density of states per or- 
bital po = 1/(2D), and take Vtk = V. The dot-lead 
coupling is then embodied in the hybridization strength 
r = 7rp|yp (with p = Afpo the total density of states, 
and JV —^ oo the number of lead orbitals). When pre- 
senting results, we use dimensionless parameters deflned 
in terms of F, viz. 



r' 



u = 



u 



U' = 



U' 



(4) 



Jh = 



JL 

r ' 



h = 



r' 



The bandwidth D is naturally taken to be the largest en- 
ergy scale in the problem, and for our NRG calculations 
in practice we take D/T =100. 

To study the models described above, the central quan- 
tities of interest are the dot Green functions G'y;cr(a;; /i) 
[^ G,j,„{t; h) = -ie{t){{d^^{t),d]^{Q)})] with associated 

spectral density Dij-criiu; h) = —-ImGij-aii^] h) {9{t) de- 
notes the unit step function) . In the absence of an applied 
magnetic field, Gy;-|-(a;; 0) = Gij-i{LL>; 0), while for any fi- 
nite h the f- and 4,-spin Green functions are naturally 
distinct. In addition to these spin-resolved quantities, 
we will later make use of their spin-summed analogs, in 
particular the spin-summed spectrum 

D,,{iu;h) = ^Y.D,j.,„{io-,h). (5) 

Moreover, in the case of the 2LM some of the physics is 



better described in terms of the symmetrized combina- 
tions of dot orbitals2i 



III. ZERO-FIELD PHYSICS AND 
LOW-ENERGY EFFECTIVE MODELS 



dea = -7^{d-la- + C?2cr), doer = —^{dlcr " ^2,7) (6) 

from which fohow the 'even-even' and 'odd-odd' Green 
functions: 



Gee;.(t^) = i [Gll-a{uj) + G22:A^) + 2Gi2;.(c^)] (7) 

Goo-A^) = \ \GixA^) + G22A^) - 2Gi2;.(c^)] . (8) 

The connection between theory and experiment is 
made via the zero-bias differential conductance, G'^^i^^. 
For the models considered above, this is obtained exactly 
from the Mcir/Wingreen approach- which gives 



2p2 /■+00 

G°iT;h) = —Go duj 



did 



NT:rDssiio;h). 



(9) 

Here f{uj) ~ (e^^"^ + 1) ^ [ks = 1) is the Fermi function, 
and Dss{^, 'i) = "^ J2i j Dij{uj; h) is the spectral density 
of the fully symmetric impurity channel [i.e. Dii{uj]h) 
for the AIM and Dee{uJ] h) for the 2LM]. The dimension- 
less prefactor Gq — sin^ 29 reflects the relative coupling 
asymmetry to the right and left leads and is maximal. 
Go = 1, for equal couplings;?^ 

As alluded to above, present methods cannot give exact 
results for the non-equilibrium situation of a finite source- 
drain bias between the leads. While recen t prog ress has 
been made in addressing this (see e.g. refs. I45ll47ll48[ ) . the 
methods used are much more computationally intensive 
and thus we make the standard approximation of neglect- 
ing the Vsd dependence of the impurity self-energy. The 
result is that 

G,{T,Vsd\h) _ 



{2eyh)Go 



A^4M + (l-A)-^^-(-) 



duj 



duj 



NTTrDssiuj;h) 
(10) 



where fi,{uj) = /(oj — /i^) with /i^ = AeV^d a-nd /i/j = 
— (1 — A)eysd- The quantity A G (0, 1) controls the par- 
titioning of the voltage split eVgd between the two leads, 
with A = 2 corresponding to a symmetric voltage drop. 
In sec. |V] we use eqn. ()10p to interpret a number of re- 
cent experimental results, where in particular we discuss 
critically the agreement between this quasi-equilibrium 
approximation and experiment. 

Finally, our numerics are obtained from the full density 
matrix (FDM)^!^ formulation of the NRG^"^ using 
the Oliveira discretization schemeSi and a generalization 
of the self-energy method of Bulla et al.^ We find it 
sufficient to keep ~ 4000 states per NRG iteration, and 
employ an NRG discretization parameter A = 3. 



To put our finite-field results in context, we consider 
briefly the zero-field physics of the two models; starting 
with the AIM, which at zero field is well understood by 
a range of complementary techniques (see e.g. ref. 'l'). 

The AIM exhibits local Fermi liquid behavior for all 
r > 0, as reflected in the RG description by a single 
stable fixed point (FP): the strong-coupling (SC) fixed 



point 



30,31 



For fixed CZ/F, the dot occupancy ni = (ni) 



increases continuously with decreasing ei, starting close 
to ni ~ when ei/F 3> 1, and tending to a maximum of 2 
when (ei + C/)/r < 1. At the point ei = -C//2 the model 
is invariant under a particle-hole (p-h) transformationi 
and hence ni — 1 precisely. 

When charge fluctuations are suppressed by a large 
U /T ^ 1, the dot occupancy tends toward integer val- 
ues, increasing more-or-less stepwise as ei is decreased 
(under a gate voltage in practice, ei ex V^ato)- In the 
singly-occupied regime (rii ~ 1), the AIM reduces un- 
der a Schrieffer- Wolff transformation^"^ of eqns. (|ll3p to a 
low-energy effective Kondo model: defining a p-h asym- 
metry parameter 77 = (l-|-2ei/C/), for a fixed — 1 < 77 < 1 
and U = U/r » 1 this yields 



Hk^ 



k,cr 



k""ka 



Js-SiO) + Kj2fLfo. (11) 



where s is a spin-^ operator describing the dot spin, 
S'"(0) = Y.a,a' fL'^aa'foa' i^ the couductiou baud/lcad 
spin density at the dot, and /oo- = jfJ2k "-ko-- The Kondo 
exchange coupling J and potential scattering strength K 
are given in terms of the original model parameters by 



PaJ = 



8 



1 



nUl-V 



poK = 



V 



nUl-V 



(12) 



such that at p-h symmetry (r/ ~ 0) the potential scatter- 
ing K = 0. Away from p-h symmetry potential scatter- 
ing is non-vanishing but, from eqn. (|12p K/J — 77/4 and 
hence for fixed 77 the model is characterized by a single 
dimensionless parameter poJ. This in turn means that 
for a given 77, all AIMs in the strongly interacting U ^ I 
regime map onto the same low-energy effective Hamil- 
tonian; and thus exhibit universal scaling of their phys- 
ical properties in terms of the low-energy Kondo scale 

rK'--Dexp[-l/(poJ)].' 

The physics of the 2LM is naturally more complicated. 
We refer the reader to ref. [2J for detailed discussion, and 
merely summarize the key points here. When considering 
the model as a function of the level energies ei and €2, it 
is more convenient to work with 

x = ei + ^U + U', y = e2 + \U + U', (13) 

since it can be shown that the model is p-h symmetric 
when a; = = y, and that its phase diagram is symmetric 
under reflection in the lines y = zLx^ That the phase 



diagram itself is non-trivial reflects the occurrence now of 
two stable FPs: the SC FP again, and the underscreened 
spin-1 (use) FP of Nozieres and Blandin.— As for the 
AIM, the SC phase is a local Fermi liquid, while the USC 
phase is a singular Fermi liquid^'*'^'^ characterized by a 
free spin-i on the dot with a In 2 residual entropy. 

Close to p-h symmetry {x = = y) the dot levels 
are each singly occupied, and in the absence of coupling 
to the lead naturally form a spin-triplet. On coupling 
to the lead, this spin-1 is reduced to an effective spin- 
i by the underscreened Kondo effect, ^"^ whence a finite 
region surrounding the p-h symmetric point belongs to 
the USC phase. On moving further away from p-h sym- 
metry [in any direction in the ix,y) plane], the model 
eventually undergoes a quantum phase transition to the 
SC phase (see Fig. 5 of ref. l24l ). The transition is of 
Kosterlitz-Thouless (KT) type^^ the Kondo scale in the 
SC phase vanishing exponentially as the boundary to the 
USC phase is approached. This holds generically except 
at points of special symmetry (specifically along the line 
y — X, where the transition becomes first-order—). 

The effective low-energy model 'deep' in the under- 
screened triplet regime can be obtained by Schrieffer- 
Wolff on the 2LM, eqns. (031), valid formally for U/T > 1 
with fixed Jh/U and U' /U. The resulting model is spin- 
1 Kondo with potential scattering^ of the same form 
as eqn. (jlip but with s now a spin-1 operator and J,K 
given by 



PQj{x,y) = 



poK{x,y) = 




7r(f7+ij) 



where (with z — x ot y) the asymmetry is 



viz) 



2z 



U 



jJh 



(14a) 



(14b) 



(15) 



The characteristic Kondo scale deep in the USC phase 
D exp[—l/{pQj)] (i.e. has the same expo- 



.,56 



T, 



s=i 



nential dependence on pqJ as for the spin- 2 case). 

In direct analogy to the AIM, the ratio of iiT to J is 
a function solely of the asymmetries, conveniently ex- 
pressed in terms of a quantity ri(x,y): 



ri{x,y) 



2K{x,y) 
J{x,y) 



r]{x) +T]{y) 
_ l'n(x)-v{y)]^ 

l-r){x)-q(y) 



(16) 



Sufficiently deep in the USC phase, one thus expects 
physical properties of the 2LM to be universal in T/T-^^^ 
for fixed rj{x,y); as considered further in sec. IIVI 



IV. FIELD DEPENDENT STATICS AND 
DYNAMICS 

We now turn to our main focus: the effect of a applied 
magnetic field on the AIM and 2LM. While much is al- 
ready known for the AIM, certain aspects of its dynamics 
in a magnetic field--""— have not been fully understood, 
and in sec. IIVBI we present NRG results to clarify the 
situation. The 2LM model has been less widely studied, 
and we consider it in somewhat more detail. 



A. Magnetization 

It is first instructive to consider the magnetization for 
level i, here defined by 



mi{h) = (fiif) - (hii) 



duj [Ai;t('^) - Ai;4,(w)] 



(17) 
(18) 



This can be determined accurately using the FDM- 
NRG )^^i^^'^° the complete Fock space approach circum- 
venting known problems arising in the original NRGi^i^ 

Fig. [T] shows the total dot magnetization m{h) = 
mi{h) vs \og{h/TK) for the symmetric AIM with U — 50 
and e = -^U^ The accuracy of the FDM-NRG is 
confirmed by the clear agreement with the exact result 
known from the Bethe ansatz^^ for the Kondo model. At 
a field h ^ Tk the magnetization rises rapidly from its 
zero field (Kondo-screened) value rn(0) = 0, before turn- 
ing over to a slow asymptotic approach to saturation of 
form m{h) - 1 - [2 \n{h/T]i)]-'^ . The inset to fig.Hgives 
results for U = 30, 20 and 10, showing the inevitable de- 
viation from the universal Kondo magnetization curve at 
sufficiently high fields h > 0{T). 

We now compare this behavior to that of the two-level 
model of eqn. ([2]). The basic physics now reflects the 
destruction of the quantum phase transition occurring 
for h = 0, and its replacement by a smooth crossover. 
In terms of FPs, the spin symmetry breaking associated 
with the magnetic field renders the USC fixed point un- 
stable for all ft. ^ 0, and so ultimately all NRG flows 
tend toward a SC fixed point (now supplemented by spin- 
dependent potential scattering). 

Fig.^upper shows the total dot magnetization ?Ti(/i) = 
mi{h) + ma (ft) for U = 2U' = 4Jh == 20 and with fixed 
ei — ~\U—U' (i.e. a; = 0, see eqn. dT^ ). upon increasing 
£2 (or y) from its p-h symmetric value. In the following 
it is useful to bear in mind that on increasing 62 at zero 
fields the model undergoes the quantum phase transition 
from USC to SC at a critical e2c — -9.94. 

Curves (a) to (c) in fig. [5] correspond to £2 < e2c and 
hence the USC phase at ft = 0. At finite field these curves 
show m{h) — >■ 1 as ft — ^ 0: an infinitesimal field fully 
polarizes the free spin-^ local moment associated with 
the USC fixed point (for ft = identically, by contrast. 
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FIG. 1: Magnetization of the AIM vs hjT^. Mam: U = 
— 2e = 50. Comparison between FDM-NRG results (crosses) 
calculated via eqn. p7|l , and the Bethe Ansatz resulli^ (line) 
for the Kondo Model. Inset: U = -21 = 30 (solid), 20 (long 
dash) and 10 (short dash) corresponding to Tk/F ~ 1.8 x 
10^", 8.0 X 10"^ and 3.1 x 10"^. Deviations from universality 
occur for h > C(F). 




FIG. 2: Upper: Total magnetization m{h) = m\{h) + m2{h) 
of the 2LM with U ^ 2U' = 4Jh = 20, ?i = -^U-U' = -20 
and £2 = (a) -20, (b) -12, (c) -10.5, (d) -9.9, (e) -9.8, 
(f) —9.5, (g) —8 and (h) —2. At zero field (a)-(c) correspond 
to the use phase; here an infinitesimal field polarizes the 
underscreened impurity moment. For larger £2 the system is 
in the SC phase at zero field, and m{h — >■ 0) — >■ 0. Lower: 
Magnetization near the crossover for h — 10"^" (solid), 10"* 
(long dashed), 10"^ (short dashed) and 10~^ (dotted). The 
h — transition is marked by a dotted vertical line at e2c ~ 
—9.94... For /i — >■ 0, m jumps discontinuously at e2c while for 
finite fields this step is smeared. 



the magnetization vanishes by symmetry). On increasing 
h, the magnetization increases monotonically, crossing 
over towards m{h) = 2 on the scale h ~ T^^^ as the 
field destroys the underscreened Kondo effect and singles 
out the Sz — +1 component of the dot triplet state»^ 
As shown in ref. |24| . T^^^ increases upon moving away 
from the center of the USC phase, hence the higher field 
required to destroy the underscreened Kondo effect for 
(c) compared to (a). 

For larger level separations, £2 > e2c (curves (d) to 
(h)), the zero-field phase is SC. Here the low-field behav- 
ior more closely resembles that of fig. [TJ At zero- field the 
dot is fully screened by the lead, and remains essentially 
so until h on the order of the SC phase Kondo scale, Tk; 
above which the spin-^ Kondo effect is progressively de- 
stroyed, and m{h) crosses over to m{h) ^ 1 associated 
with a spin-polarized spin-i on the dot. As in curves 
(a)~(c), increasing the field further then causes a second 
marked increase in m{h) when the Sz — +1 component 
of the two-electron triplet state is favored. 

Fig. [5] lower shows the magnetization at various fixed 
values of h/T as a function of £2 (focussing on the vicin- 
ity of the zero-field transition at e2c)- At any finite field, 
the magnetization decreases monotonically with increas- 
ing £2, and as h/T — )■ the curve approaches the step 
function m{e2) —> 0{e2c — £2)-— In the absence of the 
field, however, the magnetization naturally vanishes, and 
hence the limit of h — > 0+ and h = are quite distinct. 

As for the spin-i Kondo effect in fig. [U the magne- 
tization deep in the USC phase (where the low-energy 
effective model is spin-1 Kondo) is a universal function 
of h/T^^^. Fig. [3] illustrates scaling of the magnetization 
for U ^2U' = 4Jh = 30, 20 and 15 at the p-h symmet- 
ric point {x,y) — (0,0) where (hi) = 1 for i G {l,2,e,o} 
(see eqn. ^). The main figure shows m(h), the total im- 
purity magnetization. Results for different values of the 
bare parameters clearly display scaling, onto a different 
universal form than for the AIMj-M^ 

Fig. [3] (inset) shows also the magnetization of the even 
and odd impurity orbitals (such that m{h) = me{h) + 
mo{h)), for the U = 20 case. The o-orbital is clearly 
polarized to a greater extent than the e-orbital by an 
infinitesimal field, and mo{h) reaches saturation more 
quickly than me{h) (or indeed mi{h)). This reflects the 
fact that the o-orbital does not couple directly to the con- 
duction band,^** only interacting with it via the e-orbital 
which couples directly; the o-orbital as such contributing 
more to the local moment than the e-orbital. 

The situation deep in the USC phase, but away from 
p-h symmetry, is illustrated in fig. 21 As mentioned in 
sec. mil universal behavior of m{h) is expected for sys- 
tems with different bare parameters, at least for fixed 
asymmetry rj{x, y) (i.e. from eqn. (jlSp the same ratio of 
potential scattering K to Kondo coupling J). To this 
end consider first the line y — —x, for all points on which 
ri{x,y) = (eqns. (|16ll5p ). Fig. S] (inset) shows the uni- 
versal m{h) at the p-h symmetric point a; = = — y 
(line) considered also in fig. [31 compared to that obtained 
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FIG. 3: Magnetization m{h/T§^ ^) for the two-level model at 
p-h symmetry (ei — e2 = —^ — U'), deep in the USC phase. 
U = 2U' = 4Jh = 30 (solid), 20 (long dash) and 15 (short 
dash), corresponding to rj|=Vr = 1.0 x 10"*^, 7.7 x 10"^ and 
6.8 X 10"*. As with the AIM, m(/i/r|=^) exhibits universal 
scaling . ' Inset: Comparison of mx{h), x £ {l,e,o} (solid, 
long dash, short dash), for the f7 = 20 case. 




FIG. 4: Magnetization away from p-h symmetry, deep in the 
USC phase. Main: m{h) — mi (ft) -I- m2{h) for the same 
systems as fig. [3] but with ri{x) — and Tj{y) = 0.5. Clear 
universality is observed. Inset: The U — 30 case on the 
y — —X line, for x = (line) and x/T = —6 (crosses). 
The resultant universal m{h) is found to be the same in both 
cases, i.e. to be independent of asymmetry. 



some distance away from p-h symmetry at x = 6r = —y 
(crosses). The two scahng curves clearly coincide. 

The main panel of fig.|3]illustrates universality for non- 
vanishing asymmetry, showing m{h) vs h/T^^^ for the 
same U — 2U' — 4Jh values as fig. |3l but now with 
r]{x) — and ri{y) — 0.5 (i.e. r]{x,y) ~ 0.29). The three 
curves scale perfectly in the universal regime, beginning 
to deviate only at high fields h ~ 0{T). Moreover, the 
resultant universal 7n(h) is found numerically to be iden- 
tical to that arising for r]{x, y) = 0, and as such thus 
appears to be independent of asymmetry rj{x, y); a result 
we have further confirmed for a wide range of ry- values. 



B. Field Dependent Dynamics 

We turn now to the field dependence of single-particle 
dynamics for the AIM and 2LM. Much is already 
known-i^iS about the former case, but it serves as a use- 
ful comparison to the 2LM and both are experimentally 
relevant. The spin- resolved impurity spectrum is first 
considered, with a two-fold focus: the field-induced re- 
distribution of weight in the Hubbard satellites, and the 
shift of the spectral maximum from zero. 

Fig.[5]shows results for Dss;i{ijj; h) and a range of fields 
/i > 0, for both the AIM (inset) and the 2LM. The level 
energies and interaction strengths have been chosen so 
that both models are deep in the Kondo regime (for the 
AIM) or underscreened triplet (2LM), and are p-h sym- 
metric (such that Dss;-\-ii-^',h) = Dss-i{—i^]h)). In each 
case the familiar three peak structure is evident: upper 
and lower Hubbard satellites due to local charge excita- 
tions on the impurity, and a central low-energy Kondo 
resonance. We denote the half- width at half-maximum 
of the h = Kondo resonance by ujk'- the low-energy 
Kondo scale, proportional to the Kondo temperature Tk- 

In both cases, increasing the applied field causes spec- 
tral weight to be redistributed from the lower to the up- 
per Hubbard satellite, corresponding to the destabilisa- 
tion of 4--spin electrons on the dot. The striking differ- 
ence between the two is that for the 2LM (main figure), 
a significant redistribution occurs upon introducing an 
infinitesimal field (e.g. h/uj^^ = 1 x 10^^), whereas for 
the AIM this occurs only when H/ujk ~ C'(l) (inset). 
This reflects directly the behavior of the magnetization 
in fig. [5] (see eqn. ^TE\\ ): the free spin associated with the 
USC FP is fully polarized by an infinitesimal field, while 
a finite h r^ ujk is required to disrupt the Kondo singlet 
associated with the SC FP of the AIM. 

The above high-frequency behavior is relatively 
straightforward compared to that at lower energies uj ^ 
ujk', as now addressed, beginning with the AIM. The 
low-frequency behavior of the AIM spectrum in a mag- 
netic field has received significant attention using vari- 
ous techniques,—"— yet there is still some disagreement 
in the literature. Here we present results from accurate 
NRG calculations, with the aim of clarifying the issue. 

At zero field the Kondo resonance at p-h symmetry 
is centered on cj = 0, symmetric to refiection about 
a; = 0, and satisfies the Fermi liquid pinning condition 
ttTDii.„(ll) = 0) = 1. Introduction of a finite h is well 
known to shift the resonance in Dn-uiuj^h) away from 
a; = and diminish its height. ■^■*= We define Ag. as the 
magnitude of this shift, as shown in the inset of fig. [51 

The spin-summed spectrum Dii(u)]h) {— 
i ^^ Dii;o-(w; /i)) is distinct from the individual 
_Dii.cr(w; /i), since the a =t and | Kondo resonances 
are shifted in opposite directions by the field. At 
sufficiently high fields, the two resonances are far apart 
and Z?ii(w) contains two peaks separated by 2A (see 
fig. [HI inset). As h is reduced, these peaks approach 
each other and are known^i to coalesce at a field we 
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FIG. 5: Redistribution of spectral weight in the Hubbard 
bands. Main: 2LM with U = 2U' = 4Jh = 20 at p-h sym- 
metry, with hjuj^^ = (dotted), 1 x 10~^ (dashed) and 1 
(solid), Wx ^ being the Kondo scale defined as the HWHM 
of the zero-field Kondo resonance. Inset: AIM with U — 20 
also at p-h symmetry with h/ujK ~ (dotted line), I x fO~® 
(dashed) and I (solid). The result for h/ujK = 1 x 10~^ is 
coincident with that for zero field on the scale shown. 
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FIG. 7: (Color online) Universality in the field-dependence of 
the spectral maximum, Ao-(/i), in Dii;o-(cj; h) for the symmet- 
ric AIM. See text for discussion. Horizontal arrow indicates 
the exact low-field asymptote Ac/h = |. Dashed line shows 
the high-field form alog{h/uiK) + c. Solid line shows A(/i), 
the position of the maxima in the spin-summed Dii(a;; h); the 
vertical arrow indicates the exact he- 
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FIG. 6: Kondo resonance splitting in the spin-summed 
Dii{ijj;h) for the AIM, on application of a magnetic field. 
U = 20, e = -10 and (top to bottom) h/hc = 0, 1/2, 1, 4/3, 
2 and 3 with he = 0.27a;if and ujk/T = 1.7 x 10"^. Inset: 
Kondo peaks in Dii{uj;h) (solid) and 1)114(0;; /i) (dashed). 



denote he (fig- El main). Our FDM-NRG calculations 
yield a universal value he — Q.27u)k in the Kondo 
regime {U /T ^1). In terms of the quasiparticle weight 
Z = [1 — 9S]ii(w = 0)/9w]^^, easily extracted from 
FDM-NRG results for the dot self-energy I]ii(w), we 
obtain he — 0.40Zr. This is in good agreement with 
the exact result of ref.lH, hc/ZT = 0.39 . . .. 

We have performed accurate NRG calculations to de- 
termine the universal scaling behavior of A.^/h and A/ft, 
as a function of /i/wk- Before discussing these results 
it is worth explaining the calculational procedure itself. 
We find that to calculate A^ and A accurately over a 
wide range of h/uJK, it is necessary to combine results 
from different values of U . For a given C/, one cannot ob- 



tain universal results for arbitrarily high h/ujK, because 
universality arises only when h is much smaller than the 
non-universal scales F and U . Since ujk is small but fi- 
nite for a given U , there will always be a (large) h/uJK at 
which h itself becomes comparable to the non-universal 
scales and the results then deviate from universality. 

Since wk decreases exponentially with increasing U , 
this might suggest working with a very large U , for then 
one can reach very high values of h/uK before h itself 
becomes non-universal. However this is subject to a sec- 
ond problem, at the opposite end of the field scale. The 
energies that enter the Hamiltonian involve combinations 
of ft, U and F, and the double-precision arithmetic used 
in NRG thus places a lower limit on the size of ft. relative 
to U and F. If w^ is too small, then low values of h/^^K 
shift the dot energy levels by so little that they cannot 
be accurately represented in double precision. 

As such, for a given JJ there is a range of fields encom- 
passing in practice around 4-5 orders of magnitude, over 
which the universal scaling curve can be determined by 
the NRG. By combining results for different values of U 
the full scaling curve can then be built up, and by choos- 
ing Us such that the calculations overlap one can obtain 
a measure of the accuracy of the calculation. 

The points in fig. [7] show the resultant A^{h)/h ob- 
tained from a series of NRG calculations for U = 20, 40, 
60 and 100. Results for different values of U indeed over- 
lap when plotted vs h/uJK, indicating universal scaling 
behavior. At low field Ac/ft — >■ 4/3 as ft — )■ 0, recover- 
ing the exact result from Fermi liquid theoryi^i^ The 
splitting Acr/ft increases with ft/cjK, undergoing a rapid 
crossover around ft/wx ^ 1 and tending asymptotically 
to the limiting form A^/ft ^ alog^h/ux) + c, which be- 
havior agrees with results obtained from the local mo- 
ment approach 1^ 
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FIG. 8: (Color online) Field-induced shift of the under- 
screened Kondo resonance in the ph-symmetric 2LM model. 
Left panel: the spectrum _Dee;|(aj) for h — (short dash line), 
and for h/u^^ = 0.366 with two different broadening pa- 
rameters: b = 0.7 (dashed line) and h — 0.2 (solid line). Right 
panel: Ao-//i as a function of /i/oj^ ^ for ^ — ^■'^ (points), 
h = 0.4 (short-dashed), h — 0.2 (long-dashed) and 6 = 0.1 
(solid), as discussed in text. The results appear to be converg- 
ing to Ao-/ft ~ 2 for all h/ijj§^^\ this limit is marked as a ver- 
tical dotted line in the left panel for the case /i/ojk ^ = 0.366. 



We believe the 1ow-/i/wk behavior of the points in 
fig- El to be numerically exact, having repeated our 
calculations significantly more accurately and obtained 
the same results. The numerics also agree with recent 
NRG calculations^ performed in the narrow region 
--0.9 < log(/i/wK) < 0.6. As the field (and hence 
location of the spectral maximum) increases further, 
however, it becomes progressively more difficult to 
obtain accurate NRG results for Ag-. This is a direct 
consequence of the broadening procedure employed 
to obtain NRG spectra: broadening is necessarily 
performed on a logarithmic scale due to the inherent 
logarithmic discretization of the technique, so sharp 
spectral features at finite frequencies become increas- 
ingly difficult to resolve as they move away from w = 0. 
The problem can be resolved to some extent by using 
z-averaging^ and calculating the self-energy directly^^, 
but presently available computing power limits the 
extent to which this approach can be pushed. In fig. [7] 
the points were obtained by averaging results from 10 zs, 
with a broadening parameter— 6 — 0.1. Increasing the 
number of zs to 20 and working with b = 0.07 and 0.05 
gives the dot-dashed and dot-dot-dashed lines in fig. [T] 
The results are clearly sensitive to the broadening at 
high fields, although in each case they show qualitatively 
similar high-field behavior. 

Fig. |S] (right) shows analogous results for the spectral 
shift in iDee;o-(w; h) for the two-level model at p-h sym- 
metry. Here we find the results to be even more sensitive 
to the NRG broadening procedure. The points show the 
splitting obtained from averaging 10 zs with b = 0.7 (us- 



ing four different bare values of U /T as before) , while the 
short-dashed, long-dashed and solid lines are from 20 zs 
with 6 = 0.4, 0.2 and 0.1, respectively. As the accuracy 
of the calculation increases the splitting appears to be 
approaching IS.„ /h = 2 for all h, in marked contrast to 
the behavior of the AIM in fig. [71 

To pursue this further, the left-hand panel of fig. [8] 
shows Dee-i(uj]h) for a representative low-field case, 
h/uj^^ ~ 0.0366: the long-dashed line shows the spec- 
trum obtained with broadening b = 0.7, while the solid 
line shows the & = 0.2 result. The figure clearly illustrates 
the sensitivity of the finite-/i spectrum to the value of b 
employed, and in line with our conjecture above it ap- 
pears that in the limit 6 — )■ the peak position would lie 
at 2h (marked as a vertical dotted line in the figure). 

The short-dashed line in fig. [8] shows also the corre- 
sponding h = spectrum for comparison, the form of 
which (including its zero-frequency cusp) has been dis- 
cussed previouslyj2Mi It is reasonable to conjecture that 
the finite- /i spectrum has a qualitatively similar form but 
shifted so that the cusp occurs at w = 2h, although at 
present is is not possible to confirm or refute this using 
currently feasible NRG calculations. 

We conclude here with a point pursued further in 
sec. El Our discussion has concerned purely the ques- 
tion: how does the equilibrium spectrum evolve with mag- 
netic field in single and two-level dots? Here we have 
deliberately not related the equilibrium spectrum to the 
finite-bias conductance, because eqn. (jlOp is approximate 
and (as shown explicitly later in relation to recent ex- 
periments) can give quantitatively wrong results for field 
strengths in excess of a few Kondo scales;^ The figures 
shown here should not therefore be translated naively 
into quantitative predictions of conductance splittings. 
The only predictions for experiment that can currently 
be made with real certainty are those involving the zero- 
bias conductance. These are now discussed. 



C. Zero-bias conductance 

Fig. [S] illustrates universality in the zero-bias conduc- 
tance for the 2LM, as functions of h/uif^^ and T/uj^^ . 
For specificity we consider the p-h symmetric point 
(x, y) — (0, 0) (which applies also to points along the line 
y — —X deep in the USC phase, see sec. lIIII ). The inter- 
actions are set a.t U = 2U' — 4Jh, and different values of 
U = U/T are considered. For fixed h/uifJ^^ (the values 0, 
0.1, 1 and 10 are shown explicitly in fig. IQ]) the zero-bias 
conductance ^^{T, h) is seen to be universal in T /uj^^ , 
as evident from clear scaling collapse of the different U 
curves. Scaling naturally breaks down at non-universal 
scales T ~ min(r,t/), and for T ^ [/ the curves show 
peaks associated with incoherent sequential tunneling. 

Notice that at finite temperature for a given, suffi- 
ciently large h/uj^^ (in excess of ~ 0.1 in fig. JH]) there is 
a universal peak in the zero-bias conductance at a tem- 
perature T r^ h. This is analogous to the peak at finite 
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parameters), reflecting the lack of adiabatic continuity of 
the use phase to the non-interacting limit. 

On applying a magnetic field, the analysis of ref. [2J 
readily generalizes to the case of broken spin symme- 



try. Now one has 2TTrDf,e;cr{0] h) 
phase shifts for a =t, I of form 



sin Sa, with separate 



Sa ~ Trrii, 



Ilo 



where 



Ir.rr = Im Tr 



dw — G^{uj), 

OOJ 



(21) 



(22) 



FIG. 9: Scaling of zero-bias conductance for the 2LM at finite 
temperature and field. Top to bottom: h/ujj^^ — 0,0.1,1 
and 10. [/ = 2(7' = 4J with U = 30, 20 and 15 (red, blue 
and green) corresponding to cjk/F — 2.06 x 10~^, 1.64 x 10~^ 
and 1.47 x 10"^ Inset: G°{T,h) vs h/Lj%=^ for (7 = 20 and 
T/ujIT^ = 0,0.1, 1, 10. G'i{T,h = 0) vs T/ujfr^ (dashed), as 
in main figure, is shown for comparison. 



frequency in the T — equilibrium spectrum, sec. IIVBI 
Yet here the peak exists in a quantity that is both di- 
rectly measurable by experiment and calculable exactly 
by theory. Until theory is able to capture accurately 
the non-equilibrium conductance as a function of source- 
drain bias, we suggest that the field-dependence of this 
peak in the zero-bias conductance, and more generally 
the h- and T-dependence of G° , be touchstones by which 
the universality of experiment is assessed. 

The inset to fig. [5] shows another way of viewing the 
universal conductance curves. Here we fix T/oj^^ (at 
values 0, 0.1, 1 and 10, top to bottom) and vary the 
magnetic field h/ujj^^ over many orders of magnitude. 
Notice that there is no incoherent peak at large h, in 
contrast to that in the T-dependence for fixed h. This is 
because at large fields the dot is completely spin polar- 
ized, and its conductance thus weak. 

Before moving to particular experiments, we consider 
specifically the T = zero-bias conductance. For h = 
this is related to the scattering phase shift, (5, via eqn. ([9]) 
and the relation 27rrZ?ee(0) = sin 5. In previous work^l 
we derived an exact Friedel-Luttinger sum rule 



^ — ^ ^iinp 



(19) 



relating 5 to the excess charge induced by the impurity, 
'^imp (equivalent to (fii -I- ^2) in the infinite bandwidth 
limit), and the Luttinger integral II defined by 



It, = Im Tr 



d^ ^^G{lo). (20) 

0(jJ 



We showed^i that while I^ — a.s usual for the screened 
Fermi liquid phase, \Il\ — 7r/2 is by contrast charac- 
teristic of the use phase (regardless of the bare model 



and such that^^ (via eqn. ^) 



Gg(T^O) 

{2eyh)Go 



~ 5 [sin^ Sf + sin^ 61] 



(23) 



As mentioned in sec. lIVAl the USC fixed point is unstable 
for any h > 0, which means that all NRG flows terminate 
at the SC fixed point. Since the SC fixed point is char- 
acteristic of adiabatic continuity to the non-interacting 
limit, for any finite h one would expect the two Luttinger 
integrals iLa to vanish. This we have indeed confirmed 
by direct numerical calculation. 

The Luttinger integrals thus change dis continuously on 
introducing an arbitrarily small magnetic field at any 
point within the USC phase. One naturally then won- 
ders whether this has consequences for the conductance. 
To answer this one can write the phase shifts in terms of 
the excess charge and magnetization^^ defined by 



^ — ^imp,t ' '^imp.J, 
'^ ^ ^imp.t '^imp.J, 



(24) 



(with ji — {hi + 77,2) and m — mi + TO2 in the infinite 
bandwidth limit), from which eqns. (J21I23P give 



Gg(r = o) 

{2eyh)Go 



= i [1 - cos{nn{h)) cos{T:m{h))] (25) 



for any point in the (x, y) plane when h > Q. But at 
points corresponding to the USC phase at zero field, 
m{h = 0+) — 1 (see e.g. fig. [S]), i.e. it too jumps dis- 
continuously on introducing an infinitesimal field. Sub- 
stituting this into eqn. ([25]) gives 



Gg(r = 0) 

(2eV/7)Go 



7rn\ 



{h = 0+) 



(26) 



which is precisely the conductance obtained— in the USC 
phase for h = 0. In other words, although both the 
Luttinger integrals and magnetization change discontin- 
uously in the USC phase on applying a field — and hence 
the cases h = and h = 0~^ are different — the conduc- 
tance itself contains no signature of these abrupt changes. 
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Configuration 



U/T 



Tj^/K 



I 
II 



8.0 
7.2 



0.020 
0.017 



0.7 
0.65 



0.2 
0.3 



TABLE I: Parameters obtained for the two dot configurations 
of ref. |4^, by comparison to NRG results. 



V. EXPERIMENTAL RESULTS 
A. Semiconductor quantum dots 

Turning now to experiment, we begin by considering 
the work of Liu et. al. in ref. 1441 where the magnetic field 
dependence of the spin- i Kondo effect was measured in a 
GaAs device. In the experiment the gates were adjusted 
to produce two different realizations of a quantum dot 
from the same device, referred to as configurations I and 
II, with different dot-lead tunnel barriers.— 

We adopt the simplest theoretical model of the device: 
the single AIM in eqn. (m, and parameterize it using 
experimental data4i At T = the equilibrium model is 
characterized by the two dimensionless parameters''^ U/T 
and ei/r, with the experimental C/ = 1.4 meVji^. The 
level energy ei is as usual taken to depend linearly on the 
applied gate voltage Vg. we write ei + U/2 = aeSVg, with 
SVg the difference between the experimental gate voltage 
and its value in the center of the Coulomb blockade valley, 
and a a dimensionless constant. Finally, comparison to 
experimental splittings at finite bias requires the dinrcn- 
sionless quantity A (sec.|lT|, that controls the partitioning 
of the source-drain bias Vsd between the leads. 

The values of U/T and a appropriate to experiment 
could in principle be obtained by comparing experimental 
and theoretical curves for Tj^/T^ versus ei + U/2 over a 
sufficiently wide SVg range, where T^ is the Kondo scale 
at the center of the Coulomb valley (i.e. ei -I- U/2 = 0). 
We find however that the range of available data in ref.|4J 
is insufHcient to determine U/T reliably in this way, since 
near the middle of the Coulomb valley where the exper- 
imental results have been obtained, the functional form 
of the theoretical Kondo scale depends only on the ra- 
tio a^U/T, and hence U/T and a cannot be separately 
obtained. We have therefore used both the zero- and 
finite-field behavior to parameterize the model, choosing 
the best values of U/T, a and A to agree with the avail- 
able experimental data. After analyzing a wide range of 
parameter space we obtain the values shown in Table ID 

Before showing our NRG results, we comment further 
on the origin of these parameters. The values of U/T and 
A were determined first, simply by optimal fitting to the 
finite-field data at the center of the Coulomb blockade 
valley (shown in fig. 1101 below), using the approximate 
eqn. (fTO|) . Then to obtain a, we compared the experi- 
mental Tk/T^ versus gate voltage to our corresponding 
theoretical results, themselves taken from explicit NRG 
calculations. We observe that the values of a so obtained 
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FIG. 10: (Golor online) Splitting of the Kondo conductance 
peak on application of a magnetic field, at the center of the 
Coulomb blockade valley (ei — —U/2). The open squares and 
filled circles are experimental data;^ while the blue (solid) 
and red (dot-dashed) lines are NRG calculations for U/T — 8, 
A = 0.7, and U/V = 7.2, A = 0.65, respectively. The dashed 
line is Ac ~ 0, and the dotted line is Ag = ksT with the 
experimental T — 55 mK. 



are in line with the experimental estimate^"* a ~ 0.024. 
Noting that the quantity denoted 'F' in ref. |4^is 2F here, 
the ratios of U/T determined therein are 5.3 and 4.0 for 
configurations I and II, respectively. Our values are a lit- 
tle larger than these, contributory factors being: a) that 
in determining F from the widths of the charging peaks 
one must bear in mind their many-body broadening j^i^ 
which typically gives them a half- width at half-maximum 
of around 1.5 — 2r (rather than F, as arises in the non- 
interacting limit); and b) fitting Kondo scales to the Hal- 
dane formula used in ref. |4J underestimates U/T, since 
it applies asymptotically in the limit U/T ^ 1. Given 
the U/T, and the experimental U — 1.4 meV''^, we then 
calculate the Kondo temperatures T^ as shown in Ta- 
ble [H with T^ defined (as in experiment^) such that 
GO(f]^,0)/Go = e'^/h at ei = -U/2. Given the sensitiv- 
ity of absolute Kondo scales to the bare model parame- 
ters, our values are in good agreement with the experi- 
mental estimates of 0.3 K and 0.63 K (configurations I 
and II respectively) 4^ 

To add further support to these parameters we note 
that a consistent, independent determination of A can be 
obtained from the experimental conductance map. Fig. 
la of ref. |4J- The slopes of the diagonal sequential tun- 
neling peaks, when plotted with Vgd as the horizontal 
axis, are readily shown^ to be proportional to 1/A and 
1/(1 — A), and hence their ratio yields A/(l — A). From 
the experimental datai^ we extract A « 0.7, in agreement 
with our values listed above. 

With the parameters thus chosen, fig . [TOl compares the 
size of the peak splittings Aq (defined as half the peak 
to peak splitting in the finite-bias conductance, using the 
notation of ref. 44) from theory - using the approximate 
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eqn. pO)) - and experiment;^ both obtained in the centre 
of the Coulomb valley (ei = —U/2) We have plotted the 
data in the form employed in ref. H, subtracting the Zee- 
man splitting q^bB ~ 2h from the actual splitting Aq 
to emphasize the deviation of the two. The open squares 
and filled circles are the experimental data for dots I and 
II respectively (as in Fig. 4 of ref. \4M . while the blue 
(solid) and red (dot-dashed) lines are the correspond- 
ing theoretical results. The latter have been obtained 
at T = 0: we note that the experimental splittings are 
generally somewhat in excess of T (see the dotted line in 
fig. nop , and hence temperature does not play an impor- 
tant role in the analysis. Note also that for fields very 
close to the coalescence point where Ac vanishes (i.e. 
when the splittings approach the dashed curve in fig. [T0|) 
it is difficult to extract the precise value of the splitting, 
and hence we show only the sections of the curves for 
which the splitting can be determined reliably. 

The agreement between theory and experiment is very 
good. In both cases theory reproduces well the low- 
field splittings, and the curves track the experimen- 
tal results up to fields of around 5T, corresponding to 
Kondo peak splittings of around 1-2T^. At higher fields 
the theoretical curves certainly deviate from experiment, 
which we take to be a sign of the breakdown of the 
quasi-equilibrium approximation in eqn. (jlOp . Recently 
Schmitt and Anders have extended their non-equilibrum 
scattering-states NRG approach to the Anderson model 
in a magnetic field^'^^ this approach offers a promising 
means of determining peak splittings out of equilibrium, 
and further work comparing its predictions with those 
of eqn. PH)) should help to establish the regimes of the 
model where non-equilibrium effects play a large role. 
One interesting question to be pursued here is the effect 
of left-right asymmetry in the coupling to the leads^ 
since within the quasi-equilibrium approximation this af- 
fects only the dimensionless Go in eqn. (TTOl) and thus 
simply rescales the conductance uniformly. 

The final point to note here is that theory reproduces 
the crossing of the two curves identified in the experi- 
ment. We find this to be entirely a consequence of the 
slightly different As for dots I and II: if one repeats the 
calculations with equal As, the curves do not cross. This 
in fact is a special case of a more general finding: for a 
given A, our calculations show that curves with different 
U/T never cross (even if one or both curves correspond 
to non- universal parameter regimes) . Hence the crossing 
of the two curves should not be taken'*'* to indicate the 
breakdown of universal scaling per se. 



B. Carbon nanotube quantum dots 

We now turn to an analysis of the experiments of Quay 
et al.^ in which the magnetic field dependences of both 
spin-i and spin-1 Kondo effects were measured in differ- 
ent Coulomb valleys of a carbon nanotube device. 



1. Spin-^ Kondo valley 

In the spin-i valley, conductance maps were 
obtained*"^ at zero and finite-S as a function of gate and 
source-drain biases, and the evolution of the finite-bias 
conductance was also measured as a function of B. The 
splitting of the Kondo resonance at finite bias was com- 
pared to various theoretical predictions in the literature, 
the level of agreement being rather poor 4^ In this section 
we explain why the experiment did not recover the ex- 
pected behavior. First, we again parameterize the model 
from zero-field experimental data. 

As before, the spin-i Kondo effect in experiment is 
captured well by the Anderson impurity model, eqn. ([l}. 
By comparing to the experimental conductance maps in 
ref. i43 we find the value U/T — 8.5 gives optimal agree- 
ment with the experimental data at both zero and finite 
fields. The value of A ~ 0.58 can separately be extracted 
as described in the previous section,— and the exper- 
imental U can be determined from the Coulomb peak 
position in Fig. 2(d) of ref. |43: it is readily seen to be 
approximately 2 meV, and hence F ~ 0.24 meV. From 
NRG calculations at T = we then find f °. == 0.094F ~ 
0.022 meV ~ 0.25 K, lower than the experimentally es- 
timated value of 2 K. This means that the temperature 
of the device (T = 352 mK)*'^ is then on the order of the 
Kondo scale, rather than being somewhat less than it. 
We believe this to be more consistent with experiment, 
as now explained. 

The magnitude of the experimental Kondo scale can 
be gauged by inspection of Fig. 2(d) of ref. [43. If these 
results were obtained at a temperature somewhat below 
T^, the Kondo resonance would hardly be eroded by tem- 
perature, and instead one would naturally attribute the 
diminution of the zero-bias conductance from the uni- 
tarity limit of 2e'^/h to the asymmetry of the left and 
right dot- lead couplings (manifest in a Go ~ 0.4). But 
it is then difficult to explain the heights of the Kondo 
resonance relative to that of the Coulomb peaks since 
(with the caveat that eqn. PU)) is approximate) we would 
expect^ the latter to be around a quarter of the height 
of the former for T ^ T^. We believe it much more 
likely that T^ is closer to the temperature of the device, 
eroding more the Kondo resonance and thus reducing its 
height to something closer to that of the Coulomb peaks. 

Moving on to our conductance results. Fig. [TlTa) 
shows the theoretical conductance map, calculated from 
eqn. ([TU)) . to be compared to Fig. 2(a) of ref.Hj. The gen- 
eral agreement is good; the theory reproduces the intense 
sequential tunneling peaks when the dot level is resonant 
with one of the lead chemical potentials, the somewhat 
weaker Coulomb diamond, and the narrower Kondo reso- 
nance at zero-bias near the centre of the Coulomb block- 
ade valley. Fig. [TTT b) shows the effect of switching on a 
magnetic field h/T — 0.5: again, the qualitative agree- 
ment with experiment is very good, including now a clear 
'ellipsoidal' ring around the centre of the Coulomb block- 
ade valley resulting from the splitting of the Kondo res- 
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FIG. 11; Spin-i conductance maps at zero and finite magnetic 
field: to compare witii fig. 2 of ref. |43. An AIM is used with 
tj = 8.5 (giving u^kJV ~ 0.15), with T = lok and Go = 1. 
(a) Zero field conductance: a clear Kondo ridge is seen at 
zero bias voltage, (b) A finite field, h — 0.5, splits this ridge 
into twin peaks away from zero bias, (c) Evolution of the 
Kondo peaks with field at the middle of the Coulomb valley 
(p-h symmetry), (d) Vertical slices through c) with (top to 
bottom) h = 0, 0.1, 0.2, ..., 0.9, offset by 0.02 2e^/i"^ per slice. 



onance by the field. 

The field dependence is shown in more detail in 
fig. fTTTc) and (d), which both show the field-dependence 
of the conductance in the center of the Coulomb block- 
ade valley. We again recover the key features of the ex- 
periment (Fig. 2(c) and (d) of ref. |43I ). For fields suf- 
ficiently small compared to the zero-field Kondo scale 
{h/T < 0.2) the Kondo resonance remains intact, while 
for larger fields it is progressively split and ultimately de- 
stroyed with increasing h, eventually leading to a region 
of almost zero conductance around zero bias. Comparing 
the slices through the data in fig. [TlT d) to those of the 
experiment, we again observe good qualitative agreement 
between the two. 

We should point out at this stage that the value of 
h/T = 0.5 chosen in fig. [TlT b) corresponds, in physical 
units, to a field of about 2 T, around half the experimen- 
tal field. This again we attribute primarily to the break- 
down of the quasi-equilibrium approximation eqn. (llOp at 
fields larger than a couple of Tk: as seen earlier in fig.flUl 
the approximation tends to overestimate the splitting at 
these high fields, and hence a smaller h/T must be used 
in the calculation to obtain the same absolute splitting as 
the experiment. Based on the comparison of the previous 
section, noting the theoretical value here of T^ = 0.094r, 
we would estimate that the quasi-equilibrium approxima- 



tion begins to break down for this experiment at fields 
B> 1 T. 

While the latter means we cannot compare quantita- 
tively our NRG predictions at finite field to those of the 
experiment over the whole range of fields measured, we 
nonetheless believe the parameterization of the experi- 
ment to be reliable at low magnetic fields. This allows 
us to make order-of-magnitude predictions that explain 
the significance of the experimental results and the rea- 
son for the apparent disagreement with theory^ as now 
explained. 

To summarize the analysis of ref. \^ first the split- 
ting of the Kondo peak with field was extracted from the 
experimental data and plotted versus B. It was found 
that half the splitting tends to the form i5 = gfi^B at 
high field (with g ~ 2.07). Direct comparison was then 
made between the full field-dependence of the splitting 
obtained from several theories, and experiment. 

We point out that there are two basic problems with 
making this comparison. First and foremost, if one is in- 
terested in the universal form of the Kondo splitting, the 
experimental parameters need to satisfy both U/T 3> 1 
and h <C min(r, U). The former condition is necessary to 
ensure that the experiment is well-described by an effec- 
tive Kondo model at low energies, and arises because the 
Schrieffer- Wolff transformation that maps the full Ander- 
son model onto the Kondo model is formally valid in the 
asymptotic limit U/T ^ 1. The latter condition defines 
what is meant here by 'low energies': even if U/T is large, 
the effective Kondo description will always break down 
at energies of the order of the non-universal scale F, and 
the results on such an energy scale will simply not show 
universal Kondo form. 

One could argue that the U/T ~ 8.5 here is suffi- 
ciently large for the experiment to be well described by 
a Kondo model at zero field, although we believe this to 
be a somewhat more borderline case. The main prob- 
lem however is that the experimental U and F are too 
small for the high-field results to be universal. This can 
in fact be seen directly from Fig. 2(d) of ref. |4J (see 
also fig. [TlTb.d)): even at moderate fields of 2-3 T the 
Kondo ('Zeeman') peaks are already overlapping signifi- 
cantly the non-universal Coulomb peaks. More formally, 
to be universal for some given /i/Tk requires h/T <C 1; 
the experimental U = 2 meV^ and F ~ 0.2 meV as 
above, whence 5/iB-B "C F when B <C 2 T. 

The second problem is that the predictions for the 
theoretical conduct ance^^^i^ have all been made us- 
ing (either explicitly or implicitly) the approximation of 
eqn. (jlOp . rather than from a full-blown non-equilibrium 
approach. Even when we use the appropriate non- 
universal parameters in our NRG calculations, the com- 
parison to both the present experiment and that of the 
previous section suggests that eqn. pU)l is quantitatively 
reliable only for fields smaller than a few Kondo scales,— 
and even then is strongly dependent on the value of A. 
Until non-equilibrium approaches such as the scattering- 
states NRG^ become more feasible, the quantitative, 
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universal form of the conductance splitting for h ^ T^ is 
an open question; one should certainly not expect a pri- 
ori to obtain quantitative agreement between experiment 
and eqn. ([TU]). 



2. Spin-1 Kondo valley 

Finally we consider the effect of magnetic field on the 
conductance of the two- level model, eqn. ©, to make 
comparison with the spin-1 Kondo valley experiments of 
ref. I43I Since there are more interactions in the two- 
level Hamiltonian than the AIM, it is obviously harder to 
parameterize the model from the available data in a fully 
systematic manner. We have thus endeavoured to choose 
physically reasonable parameter values that reproduce 
qualitatively the experimental results (c.f . those used 
for the experimental comparison of ref. |2J) ; from which 
we find U/T = 12, U' /T = 6, Jh/F = -0.5 and (63 - 
ei)/r = Ae/r = 4.5. For simplicity we take T = 0, 
A = 5 and Gq = 1. Choosing a reasonable value of 
r ~' 0.25 meV gives e.g. a charging energy [/ ~ 3 meV 
and level spacing Ae ^ 1 nieV, both of which are within 
typical experimental estimates. 

Fig. [T^a') shows the resultant splitting of the 'under- 
screened Kondo' conductance peak at a point in the USC 
phase near the zero-field USC/SC phase transition (as 
indicated by the tail of the arrow in the phase diagram 
fig. [T^c)). The figure is qualitatively similar to that for 
the spin-i Kondo peak in a magnetic field [fig.lTlTc)] but, 
as noted in the case of a stretched spin-1 molecule,-- the 
field at which the zero-bias peak is destroyed is a some- 
what smaller fraction of the Kondo scale. This naturally 
reflects the sharper USC Kondo resonance (see fig.[5]left) 
compared to the spin-i Kondo case. 
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for the chosen parameters is 



us^ ~ 0.66r ~ 0.2 mev ~ 2 K, which again appears 
roughly in line with the widths of the Kondo peaks in 
the experimental conductance maps."*^ This means it is 
perhaps misleading to refer to the basic phenomenology 
here as 'underscreened Kondo' physics, since resonance 
widths on the order of F imply the model is far from 
being well described by a effective spin-1 Kondo model 
per se. As for the semiconducting quantum dot analyzed 
previously^! it also appears that the experimental trajec- 
tory as a function of gate voltage (ei ex V^ato) just cuts 
the 'edge' of the USC phase where the USC-phase Kondo 
scale is relatively high. 

Just across the phase boundary into the SC Fermi 
liquid phase, we obtain the conductance map shown in 
fig. [T2lb). Here we have kept the interactions and Ae 
fixed, but increased ei by 3.5r from its value in fig. ll2f a') 
(the head of the arrow in fig. [T^c) gives the precise lo- 
cation relative to the phase boundary). The qualitative 
agreement between theory and experiment (Fig. 5(b) of 
ref. I43I ) is again very good. We recover all basic features 
seen in experiment)^ at zero field the conductance peaks 
at around ±2r ~ ±0.5 meV, reflecting at zero bias the 




FIG. 12: Conductance maps for the 2LM close to the h = Q 
phase transition; for U = 12, U' = 6, Jh = —0.5 and 
Ae = 4.5. To be compared to figs 5a,b of ref. i43i . 
(a) USC phase for y = -x = 2.25r {x = e^ + \U + U' and 
y = X + Ae, with ei oc Ugatc). The conductance peak splits 
on applying a field, (b) SC phase for {x,y) — (1.25r, 5.75r). 
The h — Kondo anti-resonance is 'fiUed-in' for h > 0; full 
discussion in text, (c) h — phase diagram for above param- 
eters as a function of x/T and y/T, with the arrow showing 
the 'trajectory' taken in going from (a) to (b). (d) Zero- bias 
cut through (b) (crosses), along with the total dot occupa- 
tion number n{h) (dashed), magnetization m{h) (dotted) and 
G°/{2e'^/h) as given by eqn. ^ (sohd). 



antiresonance in the equilibrium spectrum just inside the 
SC phaseJ^i^i These peaks move toward each other, cross 
and ultimately move apart with increasing field, which 
can be loosely associated with a crossing of the isolated 
dot singlet and lowest triplet states, with a finite-field 
Kondo effect taking place at the crossing point (again, 
the 'Kondo' scale here is rather large, and as such one 
cannot describe the low-energy behavior in terms of a 
pure spin-i Kondo model). We note that in our calcu- 
lations the crossing takes place at /i ~ O.TF ~ 0.2 meV 
and hence i? ^ 3 T, again in good agreement with exper- 
iment. One can also make out various weaker features in 
the conductance, parallel to the main features and again 
seen experimentally, which mirror transitions from the 
isolated dot singlet to the higher energy triplet states.— 

The zero-bias conductance is analyzed further in 
fig. fT2ld). which is a cut through fig. [UK b) at eI4d — 
(crosses are the NRG data from fig. [T2l b)). With increas- 
ing field, G'°/(2e^//i) increases from its zero-field value of 
~ 2 J passes through a maximum at h/T ^ 0.7 (as evi- 
dent from fig.[T2lb)). and decreases monotonically there- 
after. Also shown are the total dot occupation n{h) and 
magnetization m{h) (see eqn. (I^^ ). both of which in- 
crease smoothly and monotonically as the ground state 
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evolves with increasing field. At zero field the dot is in a 
mixed-valent regime, with n{h = 0) ~ 1.5 (and m — 0). 
But with increasing field the dot ground state becomes 
progressively more like the simple Sz — 1 component of 
the isolated-dot triplet, with both total charge and mag- 
netization tending to 2 for h/T 3> 1 (i.e. riimp,t — 2, 
nimp,; — 0). 

While the zero-bias conductance shown above is cal- 
culated using eqn. (|9]), and as such probes single-particle 
spectra, its field-dependence shown in fig. [T^d) should 
equally be exphcable from eqn. ([25]) (sec. lIVCp . ex- 
pressed solely in terms of the dot charge and magnetiza- 
tion. That this is indeed so is shown directly in fig. [T2\d) : 
the solid line is calculated from eqn. (|25|) , and seen to be 
in very good agreement with the direct NRG calculations. 



VI. CONCLUDING REMARKS 

In this paper we have considered in some detail the ef- 
fects of an applied magnetic field on single- and two-level 
quantum dots tunnel-coupled to a metallic lead; includ- 
ing magnetization, single-particle dynamics and conduc- 
tance, and highlighting for the two-level model in partic- 
ular the rather subtle differences between the limits h = 
and /i — >■ 0. We have used NRG to analyze critically the 
field-dependent shift of the Kondo resonances in the two 
models, providing an algorithm that in principle can gen- 
erate the universal scaling behavior over arbitrarily large 
h/ujK ranges, limited in practice only by the logarithmic 
broadening inherent to NRG. For the single-level AIM, 
calculations can now be performed sufficiently accurately 
to achieve convergence up to fields of around h ^ IOOwk; 
for the two-level model convergence is slower, but ap- 
pears to indicate a constant spectral shift A^- = 2h for 
all fields in the universal regime. 

We have also made direct comparison between NRG 
calculations and two recent sets of conductance exper- 



iments on quantum dots in a magnetic field 
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Anderson-type models for the dots to determine bare pa- 
rameters corresponding to the experimental realizations. 
Agreement between theory and experiment is found to 
be very good qualitatively (essentially all salient exper- 
imental features are captured by the models), and even 
quantitatively - see e.g. fig. [TU]- provided the system is 
not 'too far out of equilibrium'. Since NRG provides in 
essence numerically-exact results, the deviation of calcu- 
lations from experiment provides a measure of the quan- 
titative reliability of the quasi-equilibrium approximation 
in eqn. (flUl) . used throughout to calculate conductance; 
we find it typically breaks down when the field-induced 
splitting exceeds somewhat the zero-field Kondo scale. 

We have argued that neither experimental^ consid- 
ered has measured the universal conductance splitting of 
the spin-i Kondo effect; and have emphasized (sec. IV A)) 
the considerable sensitivity of the field-dependence of the 
conductance peak to the partitioning of the bias potential 
between the leads (embodied in A) - over which, to our 
knowledge, there is relatively little experimental control. 
In addition, as above, eqn. (jlOp for the conductance is ap- 
proximate out of equilibrium, and until non-equilibrium 
approaches such as e.g. the scattering states NRG"'^ are 
sufficiently developed to become the mainstay, present 
theoretical tools are limited in that respect. 

In view of the above, we suggest that more experimen- 
tal attention should be given to the equilibrium, zero-bias 
conductance. Given exactly by eqn. 0, and independent 
of A, its field dependence can be calculated exactly (see 
e.g. sec. IIVCT) . We believe it presents a better prospect 
for ascertaining universality in the magnetic field depen- 
dence of spin-^ and spin-1 Kondo effects in real quantum 
dots. 
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